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We explain the observation of clusters in the tunneling res- 
onance spectra of small metallic particles of few nanometer 
size. Each cluster of resonances is identified with one excited 
single-electron state of the metal particle, shifted as a re- 
sult of the different nonequilibrium occupancy configurations 
of the other single-electron states. Assuming the underlying 
classical dynamics of the electrons to be chaotic, we deter- 
mine the typical shift to be A/g where g is the dimensionless 
conductance of the grain. 

An interacting many-body system exhibits, in general, 
a very complicated behavior. Usually one can analyti- 
cally characterize only statistical properties of the spec- 
trum. The fact that, for high enough energies, these 
properties are very well described by random matrix the- 
ory (RMT) was first attributed to the complexity of 
the many-body system. More recently, it has become 
clear that RMT also describes single-particle quantum 
dynamics which is chaotic in the classical limit y|| . Ex- 
amples are non-interacting electrons in small disordered 
metallic grains , and in ballistic quantum dots || . Real 
systems, however, contain a large number of interacting 
particles, and a question which naturally arises is how 
does chaos in a single-particle description manifest itself 
in the properties of the true many-body problem? 

Experimental || as well as theoretical |?],|^] studies 
of this problem, have been mainly focused on two is- 
sues: the statistical properties of the ground state en- 
ergy of quantum dots as the number of electrons changes, 
and the lifetime of a quasiparticle in such structures. 
Here we consider the nonequilibrium tunneling resonance 
spectra of small metallic particles Q. These spectra 
can be measured experimentally with high precision [see 
Figs. l(a,b)] and interpreted within the Hartree-Fock ap- 
proximation. They constitute a clear demonstration of 
the interplay between many-body interactions and quan- 
tum chaos, and also provide direct information on the 
quantum chaotic nature of the system. 

The experimental system consists of a single aluminum 
particle connected to external leads via high resistance (1 
- 5 Mf2) tunnel junctions formed by oxidizing the surface 
of the particle. In Figs. l(a,b) we plot the differential 
conductance, dl/dV, of two different particles (of sizes 
roughly 2.5 and 4.5 nm) as a function of the source-drain 
bias energy eV. The spectra display three clear features: 



600 

C! 

$5 400 



a 
o 



> 

T3 



3 



> 
"O 







1 1 
- (a) 


i i i i 


d. i 









120 



200 
150 
100 
50 




2 3 4 5 
energy (meV) 



1 2 
energy (meV) 






iii 


1 1 | 1 1 ! _ 


[ " 


(c) 






1 

11 
li 


t : 
ii 

h — 






ii 

i i aI^WU i - 



0.1 



0.3 



0.2 
eV (U) 

FIG. 1. (a,b): The low temperature (30 mK) differential 
conductance dl /dV versus bias energy of small Al particles 
with volumes (a) k, 40 nm 3 (b) « 100 nm 3 (Ref. [9]). The first 
resonance is isolated while subsequent resonances are clus- 
tered in groups. The distance between nearby groups of res- 
onances is approximately the single-particle mean level spac- 
ing A. (c): Model differential conductance obtained from 
nonequilibrium detailed-balance equations: solid line - in 
the absence of inelastic processes, l/fin = 0; dashed line - 
with inelastic relaxation rate larger than the tunneling rate, 
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(1) The low energy resonances are grouped in clusters. 
The distance between nearby clusters is of order the mean 
level spacing A of the noninteracting electrons in the dot. 

(2) The first cluster contains only a single resonance. (3) 
Higher clusters consist of several resonances spaced much 
more closely than A. 

In this Letter, we explain these features as conse- 
quences of the underlying chaotic dynamics of the con- 
fined electrons. Each cluster of resonances is identified 
with one excited single-electron state, and each reso- 
nance in turn is associated with a different occupancy 
configuration of the metal particle's other single-electron 
states. The appearance of multiple resonances reflects 
the strongly noncquilibrium state of the particle. 

Our model for the system is given by the Hamiltonian: 
H = Hq + Ht + i^int- Here Ho describes the noninter- 
acting electrons in the left (L) and right (R) leads and in 
the metallic grain, 



Ho = 



E 



E 



eicjci. 



a=L,R q 



Tunneling across the barriers is described by 

H *= E E T i Q) <Q+ h.c, 

a=L,R q,l 



(i) 



(2) 



where are the tunneling matrix elements. Interac- 
tion effects are taken into account only for the electrons 
in the grain, but including screening by image charges in 
the leads. Thus 



H 



int 



5£ R 

ijkl 



(3) 



where Uijki is the matrix element of the Coulomb inter- 
action for the electrons inside the grain. We remark that 
for the small aluminum grains considered here one can 
neglect superconducting pairing since the single-particle 
mean level spacing, 1 meV, is larger than the BCS 
superconducting gap which is 0.18 mcV fic| . 

The interaction term of the electrons is generally 

approximated by Y,ijkl U iiM c i c j c k c l ~ ( e J2i 4 c i) 2 / C > 
where C is the effective capacitance of the grain. Within 
this approximation, known as the orthodox model [pi, 
the charging energy depends only on the total number 
of electrons in the dot, but not on their particular oc- 
cupancy configuration. The orthodox model is able to 
account for the Coulomb blockade jll[] , and the Coulomb 
staircase behavior of the current as the number of ex- 
tra tunneling electrons in the dot increases. It can also 
be generalized to describe features on the scale of the 
single-particle level spacing Q. However, the orthodox 
model cannot account for the clusters of resonances in 
Fig. l(a,b), since these result from fluctuations, 5U, in 
the interaction energy between pairs of electrons. Before 
discussing the origin of these fluctuations we examine 
their effect on the differential conductance of the dot. 
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FIG. 2. An illustration of transport through the metal par- 
ticle at various values of the source-drain voltage V. Filled 
single-particle levels are indicated by full circles & empty 
ones by open circles. U is the charging energy, and A is the 
single-particle mean level spacing, (a) The system at small 
voltage bias within the Coulomb blockade regime; (b) V cor- 
responding to the first resonance in Figs. l(a,b). The thin 
dashed lines indicate the energy of a level after an electron 
has tunneled into the dot; (c) V near the first cluster of res- 
onances in Figs. l(a,b). The splitting within the first cluster 
originates from the sensitivity of level i + 1 to the different 
possible occupation configurations as shown. 



We focus our attention on the (experimental) voltage 
regime where there is no more than one extra tunnel- 
ing electron in the dot. At small voltage bias, V, within 
the Coulomb-blockade regime [Fig. 2(a)], current does 
not flow through the system. Current first starts to flow 
when one state i inside the grain becomes available for 
tunneling through the left barrier, say, as illustrated in 
Fig. 2(b). As the system becomes charged with an addi- 
tional electron, the potential energy of the other electrons 
in the dot increases by U ~ e 2 /C, and some of the lower 
energy occupied electronic states are raised above the 
right lead chemical potential [in Fig. 2(b) these "ghost" 
states are shown as dashed lines]. Electrons can tunnel 
out from these states into the right lead leaving the par- 
ticle in an excited state. There is, however, only one 
configuration of the electrons which allows an electron to 
tunnel into level i from the left lead, namely all lower 
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energy levels occupied. This implies that only a single 
resonance peak appears in the differential conductance 
at the onset of the current flow through the system (bro- 
ken spin degeneracy would cause splitting of this peak) . 

The situation changes when V increases such that elec- 
trons can tunnel from the left lead into the next higher 
available state i + 1, as shown in Fig. 2(c). In this case, 
there are several possible occupancy configurations, on 
which the exact energy of level i + 1 depends. The sev- 
eral possible energies of level i + 1 lead to a cluster of 
resonances in the differential conductance of the grain. 
The scenario described above holds provided that inelas- 
tic processes are too slow to maintain equilibrium in the 
particle. 

To explicitly demonstrate the splitting of resonances 
induced by fixed fluctuations in the interaction energy 
SU, model detailed-balance equations JlJ] were solved nu- 
merically and the corresponding differential conductance 
plotted in Fig. 1(c) by the solid line. The model system 
consists of 7 equally spaced levels, occupied alternately 
by 4 or 5 electrons, in a current-carrying steady state. 
For simplicity, the tunneling rate into each level, l/r tun 
(T L ( R )(ei) in the notation of Ref. [0), is chosen to be 
uniform, and the voltage is applied by increasing the left 
chemical potential. The temperature is 1% of the mean 
level spacing A, and the variance of the fluctuations SU 
in the interaction energy is A/5. [In the absence of fluc- 
tuations (SU — 0), dl/dV consists of single resonances 
spaced by A.] 

To estimate the fluctuations in the interaction energy 
consider the Hartree term of the interaction energy, Uh ■ 
We wish to calculate the interaction energy difference 
associated with different occupation configurations of low 
energy states. Suppose that, as illustrated in Fig. 2(c), 
these differ by a single occupation number, namely, in 
one configuration the state j is empty and f is full while 
in the other j' is empty and j is full. Then 

5(7 H =ydr 1 dr 2 |^(r 1 )| 2 C/(r 1 ,r 2 )[|^(r 2 )| 2 -|^(r 2 )| 2 ] (4) 

where the index i labels an electron state other than 
j or j', U(r 1: r 2 ) is the interaction potential. Clearly 
(SUh) = 0, where (• • •) denotes ensemble or energy av- 
eraging. We are therefore interested in fluctuations of 
SUh which emerge from the non-uniform probability dis- 
tributions of the single-particle eigenstates in real space. 
To calculate (SU H ) we approximate the interaction by 
U(ri, r 2 ) ~ v^ 1 S(y-i~ r 2 ) where v is the density of states, 
then 



(6U%) = 2^ 2 J d d rd d r'C 2 (r,r'), 



(5) 



where C(r,r') = (|V(r)^(r')| 2 ) - (|^(r)| 2 )(|V'(r , )| 2 ) is the 
probability-density correlation function. For disordered 
systems it takes the form |l3| 



where a is a symmetry factor (2 for GOE systems and 
1 for GUE), is the volume of the grain, D is the 
diffusion constant, and the sum is over the diffusion 
modes <j) n (r). Introducing the dimcnsionlcss conductance 
g = hTT 2 D/L 2 A = E c /A where L is the linear size of the 
system we obtain from (||) and (||) [[u} 



(SU H )=[c- 
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where c = \Z2a^ n |n|~ 4 /7r is a constant of order unity. 
Eq. (f7j) also applies for general chaotic systems, with 
g ~ 71/ 'A where 71 is the first non- vanishing Perron- 
Frobenius eigenvalue |15|. In essence, smaller g implies 
less uniform wave functions, so fluctuations in the inter- 
action energy increase as g decreases. Experimentally we 
find g ~ 5. Unfortunately, an analytical estimate of g 
requires precise knowledge of the shape and disorder of 
the particle which we lack [[l6| . 

Within our approximation for the interaction potential 
the Fock term, SUf, is equal to —SUh, thus apparently 
SUf + SUh = 0. However, for a more realistic inter- 
action potential SUf 7^ —SUh, and moreover, the Fock 
term exists only for electrons with parallel spins. SUh is 
therefore the typical single-electron level splitting due to 
interaction. 

More generally, when M available states below the 
highest accessible energy level (including spin), are oc- 
cupied by M' < M electrons, there are (^>) differ- 
ent occupancy configurations. The typical width of a 
cluster of resonances in this case is W x l 2 c/S.j g where 
W = min(Af - M', M'). The width of a cluster of reso- 
nances therefore increases with the source-drain voltage. 
The distance between nearby peaks of the cluster, on the 
other hand, decreases as W 1 ^ 2 /("'). This behavior can 
be seen in Fig. 1(c). 

Central to our analysis is the assumption that the 
steady-state occupation configurations of the electrons 
in the dot are far from equilibrium. This condition 
holds when the rate 1/Tj n of inelastic relaxation pro- 
cesses is smaller than the tunneling rate of an electron 
into and out of the dot, l/r tM „. In the opposite limit, 
l/^m > l/rtun, the system relaxes to equilibrium be- 
tween tunneling events, and the electrons effectively oc- 
cupy only one configuration. In this case one expects each 
resonance cluster to collapse to a single peak. This be- 
havior is illustrated by the dashed line in Fig. 1(c) where 
a large inelastic relaxation rate l/ri„ = 5/r t „„ was in- 
cluded in the detailed-balance equations. 

The results shown in Fig. 1 indicate that the metal par- 
ticle in the experimental system is indeed in a strongly 
nonequilibrium state. It is useful, however, to consider 
the various relaxation processes in our system in order to 
delimit the expected nonequilibrium regime. Relaxation 
of excited Hartree-Fock states may occur due to: (1) 
electron-electron interaction in the dot beyond Hartree- 
Fock; (2) clectron-phonon interaction; (3) Auger pro- 
cesses in which an electron in the dot relaxes while an- 
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other one in the lead is excited; (4) relaxation of an elec- 
tron in the dot as another electron tunnels out to the 
lead; (5) thermalization with the leads via tunneling. The 
last two processes are small corrections since they clearly 
happen on time scales larger than the tunneling time. 

In Rcf. it was shown that excited many-body 
states of closed systems with energy e smaller than 
(g/ logg) 1 / 2 A are merely slightly perturbed Hartree- 
Fock states. In other words, the overlap between the 
true many-body state and the corresponding Hartree- 
Fock approximation is very close to unity. This justifies 
the use of our model for the low energy resonances since 
g ps 5 therefore the energy interval < e < (g/ log(;) 1 / 2 A 
contains at least the first few excited states. At high 
source-drain voltage, however, when the dot is excited 
to energy g 1 !" 1 A < e < gA, tunneling takes place into 
quasiparticle states of width e 2 /(g 2 A) Jjj. This width is 
larger than the typical separation between nearby reso- 
nances but smaller than A. Therefore, electron-electron 
scattering will obliterate the fine structure of resonances 
for high energy excitations of the dot. 

Consider now the electron-phonon interaction. The 
temperature, 30 mK, is much smaller than the mean 
level spacing, therefore, the probability of phonon ab- 
sorption is negligible, and only emission may take place. 
The sound velocity in aluminum is v s = 6420 m/sec, 
therefore the wavelength of a phonon associated with re- 
laxation of energy lu ~ A = 1 meV is approximately 
50 A, the same as the system size. In this regime, we 
estimate the phonon emission rate to be 
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where €f is the Fermi energy (11.7 eV in Al), and p is 
the ion mass density (2.7 g/cm 3 in Al). This rate is 
that of a clean metal but reduced by a factor of tA/Ti 
where r is the elastic mean free time jT?). In ballistic 
systems, r is the traversal time across the system of an 
electron at the Fermi level. Assuming ballistic motion 
this factor is of order 10 -3 . The resulting relaxation rate 
for uj = A is therefore of order l/r e _p/j ss 10 s sec -1 
which is similar to the tunneling rate 1/T tun « 6 • 10 s 
sec -1 (corresponding to a current of 10~ 10 A through the 
particle). Thus, by increasing the resistance of the tunnel 
junctions one should be able to cross over to the near- 
equilibrium regime shown by the dashed line in Fig. 1(c). 

Relaxation due to Auger process is estimated to be 
negligible. Two factors reduce this rate considerably: (1) 
it is exponentially small in w/x where w is the width of 
the tunnel junction and \ is the screening length; (2) 
interaction between electrons on both sides of the tunnel 
junction can take place only within a very limited volume. 

In conclusion, we have shown that the low-voltage 
tunneling-resonance spectrum of a small metallic grain 
reflects a nonequilibrium electron configurations each of 
which leads to a different energy of the single-electron 
level used for tunneling. Consequently, the tunneling res- 



onances appear in clusters of width A/g. Relaxation due 
to electron-phonon interaction, which becomes impor- 
tant for high resistance tunnel barriers, will collapse the 
clusters. This effect can be used to probe the electron- 
phonon relaxation rate in nanometer size metal particles. 

The research at Harvard was supported by NSF Grant 
No. DMR-92-07956, ONR Grant No. N00014-96-1-0108, 
JSEP Grant No. N00014-89-J-1023. The research at 
Cornell was supported by the Sloan Foundation, and 
was performed at the Cornell Nanofabrication Facility, 
funded by the NSF (Grant No. ECS-9319005), Cornell 
University, and industrial affiliates. 



[1 
[2 

[3. 

[4. 
[5 

ir 

[§ 
[9 
[10 

[11 

[12 

[13 

[14 
[15 



M. L. Mehta, Random Matrices, (Academic Press, New 
York, 1991) 

G. Casati, F. Valz-Gris and I. Guarnery, Lettere al Nuovo 
Cimento 28, 279 (1980). 

O. Bohigas, M. J. Giannoni, and C. Schmidt, Phys. Rev. 
Lett. 52, 1 (1984). 

K. B. Efetov, Adv. Phys. 32, 53 (1983). 

A. V. Andreev, O. Agam, B. D. Simons, and B. L. Alt- 
shuler, Phys. Rev. Lett. 76, 3947 (1996). 

U. Sivan, F. P. Milliken, K. Milkove, S. Rishton, Y. Lee, 
J. M. Hong, V. Boegli, D. Kern, and M. deFranza, Euro- 
phys. Lett. 25, 605 (1994). 

U. Sivan, Y. Imry, and A. G. Aronov, Europhys. Lett. 
28, 115 (1994). 

B. L. Altshuler, Y. Gefen, A. Kamenev, and L. S. Levitov, 



cond-mat/9609132 



D. C. Ralph, C. T. Black, and M. Tinkham, Physica B 
218, 258 (1996). 

J. von Delft, D. S. Golubev, W. Tichy, and A. D. Zaikin, 
to be published in Phys. Rev. Lett . ; R. A. Smith and 
V. Ambegaokar, :ond-mat/9609027 . 

D. V. Averin and K. K. Likharev, in Mesoscopic Phe- 
nomena in Solids, eds. B. L. Altshuler, P. A. Lee, and 
R. A. Webb (Elsevier, NY, 1991) pp. 173-271. 
D. V. Averin and A. N. Korotkov, Zh. Eksp. Teor. Fiz. 
97, 1661 (1990) [Sov. Phys. JETP 70, 937 (1990)]. 
L. P. Gor'kov and G. M. Eliashberg, Zh. Eksp. Teor. Fiz. 
48, 1407 (1965) [Sov. Phys. JETP 21, 940 (1965)]. 



Ya. M. Blanter, cond-mat/9604101 



O. Agam, B. L. Altshuler, and A. V. Andreev, Phys. Rev. 
Lett. 75, 4389 (1995). 
[16] A naive estimate of g in ballistic systems is H/rA, where 
r is the time for an electron at the Fermi energy to cross 
the system. The metallic grains of the experiment, how- 
ever, have a roughly pancake shape. Assuming diffusive 
dynamics one can show that grA/h cx {z/r) 2 where z is 
the pancake thickness and r is its radius, g is therefore 
much smaller than fi/rA. 
[17] M. Yu Reizer and A. V. Sergeyev, Zh. Eksp. Teor. Fiz. 
90, 1056 (1986) [Sov. Phys. JETP 63, 616 (1986)]. 



4 



